Dynamical Stochastic Processes of Returns in Financial Markets 
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We study the evolution of probability distribution functions of returns, from the tick data of the 
Korean treasury bond (KTB) futures and the S&P 500 stock index, which can be described by 
means of the Fokker-Planck equation. We show that the Fokker-Planck equation and the Langevin 
equation from the estimated Kramers-Moyal coefficients are estimated directly from the empirical 
data. By analyzing the statistics of the returns, we present quantitatively the deterministic and 
random influences on financial time series for both markets, for which we can give a simple physical 
interpretation. We particularly focus on the diffusion coefficient that may be significantly important 
for the creation of a portfolio. 
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I. INTRODUCTION 

In the pioneering work of Bachelier [1], it has been 
shown that the probability of price changes can be de- 
termined in terms of the Chapman-Kolmogorov equation 
and the Wiener process. Since this study was published, 
many scientific researchers in this area have suggested 
the stochastic models to explain how the probability den- 
sity function (PDF) in delay times behaves dynamically 
on rigorous statistical principles. To find the scale invari- 
ance and universality in the stochastic model is an impor- 
tant and powerful issue in the statistical analysis of tick 
data of various financial markets [2,3]. Each financial 
market has different characteristic features but they are 
expected to have some similar and scale-invariant prop- 
erties in common. We wish to focus on the findings of 
unique trends and on the universality in them. In this 
study, we mainly analyze two different financial markets, 
the KTB and the S&P 500 stock index. The former is 
supposed to open for a short term of six months while 
the latter continues without breaks except for holidays. 
It is furthermore shown that there exists a similarity in 
the correlation of returns for both of these markets. 

Several studies on the dynamics of financial markets 
have recently focused on the understanding of the fat- 
tailed PDF of the price changes. It has been shown from 
previous works [4, 5] that price changes in the empirical 
dollar-mark exchange rates is regarded as a stochastic 
Markov process for different delay times. Ghashghaie et 
al. [6] have discussed an analogy between the dynami- 
cal high-frequency behavior of the foreign exchange mar- 
ket and the short-time dynamics of hydrodynamic tur- 
bulence, and they are shown to be similar to the energy 
cascade from large and small time scales in the theoret- 
ical model of turbulence [7,8]. Furthermore, Ivanova et 
al. [9] showed the validity of the Fokker-Planck approach 



by testing the various financial markets, i.e., NASDAQ, 
NIKKEI 225, and JPY/DEM and DEM/USD currency 
exchange rates. 

Motivated by the work of Fricdrich and Ivanova 
[4,5,9], we apply the Fokker-Planck equation approach 
to two sets of financial data, in this case the KTB and the 
S&P 500 stock index. In this paper, we derive the func- 
tional form of the Fokker-Planck equation from the esti- 
mated Kramers-Moyal coefficients and present an equiva- 
lent Langevin equation for each. In Section 2, we present 
the numerical method for estimating Kramers-Moyal co- 
efficients from price returns of the real financial time se- 
ries. We present some results obtained by the numerical 
simulations in Section 3, and our concluding remarks are 
in the final section. 



II. THE FOKKER-PLANCK EQUATION AND 
THE LANGEVIN EQUATION 

In our study the return r(t), given by r(t) = 
\n[p(t + At)/p(t)] for the price p(t) at time t, is as- 
sumed to follow a Markov process. This assumption 
becomes plausible by the validity of the Chapman- 
Kolmogorov equation. To make a statistical analysis 
of financial markets, we consider a iV-point joint PDF, 
p(rN, Atpj; rjy_i, Atpj^i; r\, Ati). However, in the 
case of a Markov process, it is sufficient to consider the 
joint PDF generated by a product of the conditional 
probabilities in the following way: 



p(r N ,At N ; ...;n,A*i) 

= p(tn, M N \r N -i, Aijv-i) 



■p(r 2 ,At 2 \r 1 ,Ati) 
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FIG. 1: Plot of the time series and returns in the KTB503 
((a)and (c)) and the S&P 500 stock index ((b)and (d)). 



where the conditional probability p(x%, At 2 \xx, Ah) is 
defined by 

p(r 2 , At 2 \n, Ah) = p(r 2 , At 2 ;n, Ah)/p(r u Ah) (2) 

and evaluated for the time delay At 2 < Ah directly from 
the tick data set. In addition we consider the Chapman- 
Kolmogorov equation as a necessary condition 

p(r 2 ,At 2 \r 1 ,Ah) = J dnp{r 2 , At 2 \n, AU)p(n, Au\n, Ah) 

(3) 

which holds for any time lag Ah, with At 2 < Ah < At\. 

In order to derive an evolution for the change of distri- 
bution functions across the time scales, we should show 
the validity of the Chapman-Kolmogorov equation. From 
the validity of the Chapman-Kolmogorov equation, we 
can obtain an effective Fokker-Planck equation. This 
equation is represented in terms of the differential form 
as follows: 



^(r.At^-^iMr.At) 



-^D 2 (r,At)]p(r,At), 



(4) 

where the drift coefficient Di(r, At) and diffusion coeffi- 
cient D 2 (r, At) can be estimated directly from moments 
Mk of the conditional probability. The relations between 
coefficients and moments in Kramers-Moyal expansion 
are described as follows: 

D k (r,At) = i lim M k , (5) 



1 



M fe = — / dr (r - rfp(r ,t+ At\r,t) 



(6) 



Therefore we can derive a Lagevin equation from Eq. (4). 
d 



dr 



r(r) = Di(r(r), r) + y/D 2 (r(T),r)f(r), (7) 
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FIG. 2: Plot of the PDFs of returns for time lags At = 
1,2,4,16,32 min and day (from bottom to top) in (a) the 
KTB503 and (b) the S&P 500 stock index. The top figures 
(circles) of (a) and (b) approximately approach to the Gaus- 
sian PDF. 



where /(r) is a fluctuating (^-correlated white noise with 
a Gaussian distribution holding for (/(t)/(t )) = 2S(t — 
r). 

III. NUMERICAL CALCULATIONS 

We analyze two different databases, the KTB 503 and 
the S&P 500 stock index. The KTB 503 is a minutely 
tick-by-tick data set consisting of 13282 quotes taken 
from 1st October 2004 to 31st March 2005, and the S&P 
500 stock index is a daily data set consisting of 14027 
quotes taken from 1950 to 2005. From both data sets 
we measure the price returns r(t) for analysis, which is 
very useful to see the gain or loss of an asset in a given 
time period. Both time series are plotted together with 
the returns of the time delay At max = 32 minutes and 
32 days for the KTB 503 and S&P 500 stock index, re- 
spectively in Fig. 1. We consider a logarithmic time 



3 



x10" 3 

5 1 1 1 1 1 1 1 



1 -3 


-2 


-1 

AX, 


1 


2 


3 < 
x10" 3 


(b) 





































-0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 

AX 1 

FIG. 3: Contour plot of the joint PDF p(r 2 , At 2 ;n, Ati) for 
the simultaneous occurrence of price returns with time scales 
At 2 — 32 (min, day) and Ati = 16 (min, day) in (a) the 
KTB503 and (b) the S&P 500 stock index, respectively. 




FIG. 4: Contour plot of the conditional PDF 
p(r 2 ,At 2 \n,Ati) for S&P 500 stock index with Ati = 16 
and At2 — 32 days, (a), (b), and (c) are corresponding cuts 
for ri = -0.05, 0, and 0.05 



scale r = \n(At max / At) for the sake of convenience in 
the analysis. 

We show the empirical PDF for various time delays 
in Fig. 2 from which we can confirm the tendency for 
the PDF to approach a Gaussian PDF as the time lag 
At increases. In a comparison with the PDF of Ref. 
[4], it is shown that the numerical iteration of the effec- 
tive Fokker-Planck equation fits well the empirical PDF. 
In this paper, we restrict ourself to obtain an effective 
Fokker-Planck equation and to confirm that the approach 
presented is also well applied to the other financial mar- 
kets different from the foreign exchange market. For the 
following analysis, we calculate conditional probabilities 
by using the contour plot of the joint PDF for the si- 
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FIG. 5: Matrix representations of the conditional PDFs of 
S&P 500 stock index scaled in a jet colorbar. (a), (b), and (c) 
correspond to p(r 2 , At 2 \ri, Ati) with each (Ati, At 2 ) set to 
(8, 16), (16, 32), and (8, 32) in days in sequence, (d) indicates 
the different conditional PDF between the r.h.s and the l.h.s. 
of the Chapman-Kolmogorov equation, Eq. (3). 

multaneous occurrence of price returns, as shown in Fig. 
3. 

We can obtain approximately conditional probabilities 
by counting the simultaneous events belonging to each 
bin by which the corresponding events can be divided 
(see Fig. 4). 

Fig. 5 shows the matrix representation of conditional 
probabilities, which is a good indicator for proving the 
validity of the Chapman-Kolmogorov equation for the 
empirical data. The matrix is generated from a histogram 
with 100 discretionary steps over the range of r(t). There 
are visible deviations in the outer region of returns, which 
are probably due to the finite resolution of the statistics 
and to the small number of data set, where At max — 32 
days and 32 minutes for S&P500 stock index and the 
KTB503, respectively. 

To calculate Mk approximately we apply the condi- 
tional PDF obtained from the procedure, as shown in 
Figs. 3 and 4. We set At = In 2 in the estimation of 
Dk, and the functional r dependence of D\ and D2 is 
shown in Fig. 6. We confirmed that there is no signifi- 
cant difference in the estimated coefficients when varying 
tAt m ax, provided that the At is fixed. For all scales t and 
At, the functional dependence of D\ and D2 behaves in 
the same manner. D\ shows a linear dependence on re- 
turns, while D2 is approximated by a polynomial of two 
degrees in the logarithmic returns. As specified in Ref. 
[4] we add an exponential term to correct the delicate 
differences in the minimum values of D2 although this is 
too small to be considered significant. By analyzing the 
plots in Fig. 6, we obtained the following approximations 
for the KTB503 and S&P 500 stock index: For S&P 500 
stock index 

£>i = -0.70r (8) 

and 

L>9 = 0.48(r - 5.06 x 10~ 3 ) 2 + 7.63 x 10~ 4 cxp(-l.Olr). 

(9) 
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FIG. 6: Drift and diffusion coefficients estimated from the 
conditional PDFs. The solid curves present linear and 
quadratic fits in the KTB503 ((a) and (b)) and S&P 500 stock 
index ((c) and (d)). 

For the KTB503 

Di = -0.78r (10) 

and 

D 2 = 0.58(r - 8.21 x 1CT 5 ) 2 + 1.31 x 1CT 7 exp(-4.46r). 

(11) 

In Eqs. (9) and (11), we obtained the exponential 
term by solving the simultaneous equations from a pair of 
two conditional conditions, in this case pfa, At2\r\, Ati) 
where At\ and At 2 are chosen to satisfy the condition, 
t = In 2. The analogy between turbulence and the for- 
eign exchange market was well described and argued in 
Ref. [4]. Here we extend the argument to the stock and 
futures exchange markets. Although, as mentioned pre- 
viously, the KTB503 has a short-term expiration differ- 
ent from the S&P500 stock index, they both show good 
agreement in the form of the Kramers-Moyal coefficients. 
We here stress again that the Fokker-Planck approach 
can be a suitable method to analyze the financial tick 
data and in understanding the statistics of their returns. 
Furthermore, we present a nonlinear Langevin equation, 



FIG. 7: Plot of Z^open circle) vs ^(diamond), (a) and 
(b) correspond to the KTB503 and S&P 500 stock index, 
respectively. 

here an underlying stochastic process, by showing that 
the Kramers-Moyal coefficients vanish for k > 3 (see Fig. 
7). It can be particularly confirmed that the third mo- 
ment of Kramers-Moyal coefficients vanishes. In Eq. (9) 
and (11), the exponential terms are too small to be negli- 
gible when compared to other terms. Hence, from the es- 
timated results we can obtain a linear stochastic equation 
with multiplicative noise, as a quadratic noise Ornstein- 
Uhlcnbeck process, as follows: For S&P 500 stock index, 

-^-r(r) = -0.70r(r) + Va48r(r)/(r). (12) 
For the KTB503, 

-^-r(r) = -0.78r(r) + V058r (r)/(r). (13) 
dr 

From Eqs. (8)-(ll) we can sec that the KTB503 has a 
greater drift and diffusion coefficients than the S&P500 
stock index. This implies that the former is more dy- 
namic than the latter. Furthermore, the volatility of re- 
turns has considerable influence upon the deterministic 
trend in that the diffusion coefficient D 2 is somewhat 
comparable to the drift coefficient D\. 

IV. CONCLUSIONS 

We have showed that the Fokker-Planck approach can 
be a good method for the analysis of the financial mar- 
kets in a wide spectrum. Even for a short-term expired 
financial market such as the Korean treasury bond fu- 
tures we could obtain a well-fitted result describing the 
behavior of its PDF in time delays. We have also ob- 
tained the same stochastic process from both analyzed 
financial markets irrespective of their data acquisition 
frequency. Quantitative measurements of deterministic 
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and random influences for both futures and stock mar- 
kets imply that the futures market is more dynamic than 
the stock market. 

This paper has shown that the temporal correlation 
of the futures market can be described well in terms of 
Fokker-Planck equations. Particularly, the KTB503 is 
transacted for a short term of six months in the Korean 
financial market, but it is striking that its functional form 
of the Kramers-Moyal coefficients is approximately con- 
sistent with that of other assets. Moreover, the compar- 
ison of the drift and diffusion coefficients with various 
financial markets provides us with a numerical indicator 
for the creation of portfolios and risk management strate- 



gies. In the future, it will be of interest to compare our 
results with those of other national options. Thereby we 
can measure and compare the stability and the efficiency 
of various financial markets. 
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